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Abstract 

Oh! 

. We present the results of a calculation of the rate for the annihilation of an e^e 

pair into a jet containing an energetic photon and a single further jet at order aas- 
By comparing these results with the existing data on the photon + 1 jet rate from 
the ALEPH Collaboration, we make a first determination of the process independent 
5^ \ quark-to-photon fragmentation function Dq^^{z, iip) at order aus- We also compare 

our prediction for the inclusive photon energy spectrum with the recent OPAL mea- 
surement. 



The production of final state photons at large transverse momenta is one of the key 
observables studied in hadronic collisions. Data on high-p^ photon production yield valuable 
information on the gluon distribution in the proton. Moreover, if the Standard Model Higgs 
boson has a mass lying below the VTVT-threshold, one of the best detection channels at 
the CERN LHC is via its two photons decay mode A precise understanding of the 
background processes yielding two photon final states is therefore crucial for Higgs searches 
at the LHC. 

Photons produced in hadronic collisions arise essentially from two different sources: 'di- 
rect' or 'prompt' photon production via the hard partonic processes such as qg 97 and 
QQ ~^ 91 through the 'fragmentation' of a hadronic jet into a single photon carrying a 
large fraction of the jet energy. The former gives rise to perturbatively calculable short- 
distance contributions whereas the latter is primarily a long distance process which cannot 
be calculated perturbatively and is described in terms of the quark-to-photon fragmentation 
function 0. In principle, this fragmentation contribution could be suppressed to a certain 
extent by imposing isolation cuts on the photon, typically formulated by limiting the amount 
of hadronic energy allowed within an isolation cone around the photon. However, the match- 
ing of experimental isolation cones onto a theoretical calculation beyond the lowest order 
level is far from trivial 0, and it has been shown that a sizeable contribution to the 
isolated photon cross section from the quark-to-photon fragmentation process remains even 
after isolation cuts have been applied p. 

All estimates of the fragmentation contribution to the production of one or two final state 
photons in hadron collisions so far have had to rely on a model of the photon fragmenta- 
tion function. This function obeys a perturbative evolution equation with a non-perturbative 
boundary condition usually parameterized in the form of an initial distribution at some low 
starting scale hq. Several models, based either on an asymptotic solution of the evolution 
equation or on estimates of the non-perturbative boundary conditions inspired by Vector 
Meson Dominance models ||^ can be found in the literature. It is clear that the lack of a pre- 
cise knowledge on the quark-to-photon fragmentation function introduces a non- quantifiable 
systematic uncertainty into any theoretical calculation of photon production cross sections. 
A precise determination of the process independent quark-to-photon fragmentation function 
is therefore needed and the results obtained in the remainder of this letter should help to 
achieve this. 

A first measurement of the quark-to-photon fragmentation function has been performed 
by the ALEPH collaboration at CERN from an analysis of two jet events in which one of the 
jets contains a photon carrying a large fraction [z > 0.7) of the jet energy. In these 'photon' 
+ 1 jet events, a democratic approach is used to group the particles into jets and the photon 
is clustered simultaneously with the other hadrons. A comparison between this measured 
rate and the calculated rate up to 0{a) [|l3] yielded a first determination of the quark- 
to-photon fragmentation function accurate at this order. Both theoretical and experimental 



analyses used the Durham jet algorithm Note that z, the fraction of electromagnetic 
energy within the 'photon' jet, is in general different from the parameter considered in 



inclusive photon production |12|. In the collinear core of the 'photon jet' however, where the 



non-perturbative effects arise, the definitions of these two parameters coincide. 
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In the present letter, we report the results of a calculation of the 0{aas) QCD corrections 
to the 'photon' + 1 jet rate, performed using the same democratic approach. Together with 
the published data 0, this enables us to obtain a first determination of the process inde- 
pendent quark-to-photon fragmentation function at 0{aas)- The details of this calculation 
will be presented in a subsequent publication |13 . 



A common feature of all model estimates of the quark-to-photon fragmentation func- 
tion 1^, is the resummation of leading (a" log"'*'^ and subleading (a" log" fj,%) loga- 
rithms of the mass factorization scale fip to all orders in the strong coupling constant. The 
results of this resummation procedure seem to suggest that the photon fragmentation func- 
tion is 0{a/as) and that quark and gluon fragmentation into photons contribute at the 
same order. This view is widely held for the photon fragmentation as well as for the photon 
structure function. However, in the region of high z and where jets are resolved, this is not 
necessarily the correct procedure to adopt. For example, the distribution of electromagnetic 
energy within the photon jet of photon + 1 jet events, for a single quark of charge e^, can 
be written ]T(][, 



= 2D,^.{z, I^f) + ^P^^\z) log f ^) + . . . , (1) 

where . . . represents terms where the photon is perturbatively isolated (proportional to 
5(1 — z)) and other perturbative contributions which are well behaved as z 1. The 
non-perturbative fragmentation function is a solution of, 

-"''-Pi°) (.) + ..., (2) 



d log fXp 2n 

and Scut is determined by the jet clustering algorithm. In the Durham jet clustering al- 
gorithm and at large z, Scut ~ sz{l — zY/{l + z) ^ [|10| where px is the transverse 
momentum of the photon with respect to the cluster. In the usual approach, the solution 
of the evolution equation for Dg^^, is proportional to log(/i|'//io) where fio is some small 
scale and the conventional assignment of a power of l/a^ to the fragmentation function is 
clearly motivated. For 'physical' scale choices, //p ~ Scutfl^nd the perturbative contribution 
is suppressed relative to the fragmentation contribution, so that for z < 1, the fragmentation 
function is indeed more significant. However, at large z, we see that the transverse size of 
the photon jet cluster decreases such that Scut 0. The hierarchy Scut ~ fJ-p ^ /ig 

is no longer preserved and both contributions in eq. (|I]) are important. Large logarithms of 
{1 — z) are the most dominant contributions. Since, from the experimental point of view, 
we are required to observe isolated or almost isolated photons, this is precisely the region 
that we are most interested in. For this reason, we choose not to impose the conventional 
prejudice and resum the logarithms of hf a priori. Instead, we work within a fixed order 
framework, to isolate the relevant large logarithms. 

As mentioned above, potentially large logarithms of {1 — z) are present and need to be 
resummed. At present however, we have not accomplished this resummation and, to avoid 



^For inclusive processes, Scut s is more appropriate. 
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Figure 1: Parton level subprocesses contributing to the photon + 1 jet rate at 0{aas). 



unstable behaviour at large z, we choose a form of the fragmentation function that renders 
the calculation well behaved as 2; 1. This approach is identical to the procedure applied in 
the lowest order calculation and in the lowest order experimental analysis |0. It is thus 



possible to make a direct comparison of our next-to-leading order results with these earlier 
leading order results, and to test the stability of the perturbative fixed order approach. 

The photon + 1 jet rate in e^e~ annihilation at 0{aas) receives contributions from five 
parton-level subprocesses displayed in Fig. |l|: 

a) The tree level process 7* — > qqg'y, where the final state particles are clustered together 
such that a "photon jet" and one additional jet are observed in the final state. 

b) The one loop gluon correction to the 7* — > qq^ process, where the photon and one of 
the quarks are clustered together or the photon is isolated while both quarks form a 
single jet. 

c) The process 7* qqg, where one of the quarks fragments into a photon while the 
remaining partons form only a single jet. 

d) The one loop gluon correction to the 7* — >■ qq process, where one of the quarks frag- 
ments into a photon. 

e) The tree level process 7* — > qq with a generic 0{aas) counterterm present in the bare 
quark-to-photon fragmentation function. 

Although the 'photon' + 1 jet cross section is finite at 0{aas), all these contributions 
contain divergences (when the photon and/or the gluon are coUinear with one of the quarks. 
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when the gluon is soft or because the bare quark-to-photon fragmentation function contains 
infinite counter terms). These divergences need to be isolated and cancelled analytically. 
It is only after this cancellation has taken place, that it is possible to numerically evaluate 
the remaining finite contributions in a Monte Carlo program containing the experimental jet 
reconstruction algorithm, and ultimately to evaluate the 'photon' + 1 jet rate at 0{aas)- 

The assignment of a single power of a to the quark-to-photon fragmentation function (and 
0{aas) to the gluon-to-photon fragmentation function) carries through when factorizing 
the explicit singularities present in the real and virtual processes with photon emission 
into the fragmentation functions. Here, the leading perturbative counterterm in the quark- 
to-photon fragmentation function is 0{a) while the leading counterterm in the gluon-to- 
photon fragmentation function is 0{aas)- Consequently, the gluon-to-photon fragmentation 
function plays no role at the fixed order 0{aas)', since radiating the gluon contributes an 
additional factor of as, it first contributes at 0{aa1). In fact, all of the model estimates |^, |^ 
indicate that the gluon-to-photon fragmentation is orders of magnitude smaller than the 
quark-to-photon fragmentation function at large z and we ignore it throughout. 



Following the philosophy of we use the concept of a theoretical resolution criterion 
Sinin to define the singular regions. Within the singular regions, the matrix elements are 
approximated and the unresolved variables analytically integrated. At next-to-leading or- 
der, most QCD processes receive singular contributions when either one gluon is soft or two 
particles are coUinear - in other words a single particle is theoretically unresolved. However, 
the evaluation of the singular contributions associated with the process 7* — > gggj is of 
particular interest as it contains various ingredients which could directly be applied to the 
calculation of jet observables at next-to-next-to-leading order. For example, there are contri- 
butions where two of the final state partons are theoretically unresolved. The three different 
double unresolved contributions which occur in this calculation are: the triple collinear con- 
tributions, arising when the photon and the gluon are simultaneously collinear to one of the 
quarks, the soft/ collinear contributions arising when the photon is collinear to one of the 
quarks while the gluon is soft and the double single collinear contributions, resulting when 
the photon is collinear to one of the quarks while the gluon is collinear to the other. In 
addition, there are also single unresolved contributions from the one-loop 7* — * qq'y process 
and from the tree-level fragmentation process 7* —>■ qqg when two of the final state particles 
are collinear. A detailed derivation of each of these singular contributions, along with how 
the different single and double unresolved regions match onto each other, will be presented 
in 
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The combination of the unresolved contributions present in the processes shown in 
Fig. |l|(a)-(d) yields a result that still contains single and double poles in e. These pole 
terms however are proportional to the universal next-to-leading order splitting function 
T5| and a convolution of two lowest order splitting functions, (-P^g-* (8>-Pg^°-*). Hence, they are 



factorized into the next-to-leading order counterterm of the bare quark-to photon fragmen- 
tation function [|TB] present in the contribution depicted in Fig. |I](e), yielding a finite and 



factorization scale (/li?) dependent result |T^. 



Once the singularities have been removed, the remaining finite contributions can be 
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numerically evaluated. Often this is realized within the phase space slicing method, first 
introduced by [0 and further developed in [|14| for general next-to-leading order QCD pro- 



cesses. The original formulation of the phase space slicing method requires the analytical 
calculation of an approximated matrix element inside the theoretically unresolved regions, 
while the full matrix element is evaluated numerically only outside any singular region. This 
procedure turns out to be inappropriate where more than one particle is potentially unre- 
solved. Apart from generically double unresolved regions, we find areas in the four parton 
phase space which belong simultaneously to two different single unresolved regions. Those 
areas are not treated correctly within the phase space slicing procedure. A consistent treat- 
ment is possible in the hybrid subtraction method developed in |jl8[. This scheme attempts to 



combine the advantages of the phase space slicing procedure with the alternative subtraction 



methods [|T9|. The parton resolution parameter is used to isolate the poles, but, rather than 
assuming that the approximated matrix elements are exact, the difference between the full 
matrix element and its approximation is numerically evaluated in all unresolved regions. 

The numerical program finally evaluating the 'photon' -|- 1 jet rate at 0{aas) contains 
four separate contributions. These contributions are determined according to the number 
of partons in the final state and the presence or absence of the quark-to-photon fragmenta- 
tion function. For each final state, 'photon' -|- 1 jet events are selected using a democratic 
jet clustering algorithm and a given value of the jet resolution parameter ycnt- The "frag- 
mentation function" used in the numerical evaluation of the cross section is the sum of the 
/iiT-dependent quark-to-photon fragmentation function and a finite contribution from the 
unresolved photon region. The four contributions are: 



(A) 2 partons -|- photon 

This contribution corresponds to the process 7* — > qq'j with a hard photon in the 
final state. It is present at 0{a) and 0{aas) (due to the presence of a theoretically 
unresolved real or virtual gluon). 

(B) 2 partons -|- "fragmentation" 

The 7* qq®Dq^^ process is present at 0{a.) and 0{aas). In particular, it contains 
the finite terms corresponding to the g — 7 coUinear region at 0{a) and the finite parts 
associated with all double unresolved regions at 0{aas). 

(C) 3 partons -|- photon 

This contribution is only present at 0{aas) and describes the 7* qq'jg process where 
both photon and gluon are theoretically resolved. 

(D) 3 partons -|- "fragmentation" 

The 7* — > qqg ® Dg^^ process with a hard gluon in the final state, is only present at 
0{aas). 



Each of the contributions listed above depends logarithmically (in fact as log'^(ymin)) on 
the theoretical resolution parameter i/min = Smin/Q^, which is used throughout the calculation 
to divide the final state phase space into different resolved and unresolved regions. However, 
the physical 'photon' + 1 jet cross section, which is the sum of all theoretically resolved and 
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Figure 2: Sum of all 0{aas) contributions to the total photon + 1 jet rate for a single quark 
of charge such that ae^ = 27r, as(A^^ — 1)/2N = 27c using the Durham jet algorithm with 
Vent = 0.1, and integrated for 2; > 0.7. 



unresolved contributions, must of course be independent of the choice of ?/min- That this is 
indeed the case is illustrated in Fig. ^ where we show the sum of all 0{aas) contributions 
to the total rate as a function of ymin- For simplicity, we show only the 0{aas) contribution 
for a single quark of charge Cg such that ae^ = 27r, as{N'^ — 1)/2N = 27r using the Durham 
jet algorithm with ?/cut = 0.1, and integrated for z > 0.7. We see that the cross section 
approaches (within numerical errors) a constant value provided that y^i^ is chosen to be 
small enough, indicating a complete cancellation of all powers of log(?/min). This provides 
a strong check on the correctness of our results and on the consistency of our approach. 
The behavior at large values of T/min can be understood to be due to terms proportional 
to ?/min log^(ymin), which havc bccu neglected in the analytic evaluation of all theoretically 
unresolved contributions. It can clearly be seen that these terms are already unimportant 
at 2/min = 10^^, which will be used in all numerical studies below. 

We can now determine the quark-to-photon fragmentation function D{z,^f) up to this 
order from a comparison between the measured 'photon' + 1 jet rate and the results of our 
calculation. This function, which parameterizes the perturbatively incalculable long-distance 
contribution to the photon radiation off a quark has to satisfy a perturbative evolution 
equation in the factorization scale ^p. All unknown long-distance contributions can therefore 
be attributed to the behavior of D{z,fiQ), the initial value of this fragmentation function, 
which is to be fitted to the experimental data, at some low scale /xq- Indeed, the next-to- 
leading order fragmentation function can be expressed as an exact solution of the evolution 
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Figure 3: Comparison of the 'photon' + 1 jet rate at 0{a) and 0{aas) with the ALEPH 
data (?/cut = 0.06), including the fitted quark-to-photon fragmentation function. 



equation up to 0{aas) |13 



ae 




Ztt ^ ' \f^o J Ztt Ztt 



2N 



p{i) 



97 



27T 2 



(3) 



It should be re-emphasized that this solution does not take the commonly implemented 
(e.g. 0, H) resummations of log fip into account. As explained above, this appears the only 
consistent procedure since powers of log(l — z) may be of equal importance in the large z 
region |^ , which is the main region of interest in the present study. Finally, this procedure 
is identical to the procedure applied in earlier theoretical |jlO[ and experimental [^] analyses 
at lowest order, such that a direct and self-consistent quantification of the next-to-leading 
order effects in the 'photon' + 1 jet rate becomes feasible. 

Furthermore, since this solution is exact at the order under consideration, the factoriza- 
tion scale dependence of the 'photon' + 1 jet rate is eliminated, i.e. our results are indepen- 
dent of the choice of fip- 



A three parameter fit to the ALEPH 'photon' + 1 jet data for the z distribution. 
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Figure 4: The measured quark-to-photon fragmentation function at leading and next-to- 
leading (this study) order for a quark of unit charge. The factorization scale is taken to be 
fJ^F = Mz. 



for the jet resolution parameter ?/cut = 0.06 and as{M^) = 0.1240 yields, 



D 



NLO, 



2n 



(-n?n^)log(l 



+ 20.8(1 -z)- 11.07 



(4) 



where /io = 0.64 GeV. For reference, the lowest order fit obtained by the ALEPH Collabo- 
ration P is, 

D^«(z,/zo) = ^ (-i^;n^)log(l - zr - 13.26) , (5) 

with /io = 0.14 GeV. In both cases, the logarithmic term proportional to P^^\z) is introduced 
to ensure that the predicted z distribution is well behaved as 2 — > 1 [|l^. It should be kept 
in mind that both parameterizations are only fitted to data with z > 0.7, which is thus the 
natural limit of their validity. 

The ALEPH data and the results of the 0{aas) calculation using the fitted next-to- 
leading order fragmentation function are shown in Fig. ^ (For reference, we have also 
displayed the 0{a) prediction |1^ using the ALEPH fit of the lowest order fragmentation 
function |^.) The fit also provides a good description of the z distribution for other values 
of Vcut ||13||. This is to be expected since the non-perturbative effects are largely in the 



^At this order, all choices of the strong coupling constant are equally valid. This value of as is chosen so 
that the observed hadronic cross section is reproduced by the 0{as) calculation. 
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Figure 5: The integrated photon + 1 jet rate above z = 0.95 as function of ?/cut, compared 
with the 0{a) and 0{aas) order calculations including the appropriate quark-to-photon 
fragmentation functions. 



collinear core of the 'photon' jet which is fully contained within the jet definition. Moreover, 
it turns out that the next-to- leading order corrections are moderate for all values of ?/cut, 
which clearly demonstrates the perturbative stability of our fixed order approach. 

The resulting next-to- leading order (MS) quark-to-photon fragmentation function (for a 
quark of unit charge) at a factorization scale fip = Mz is shown in Fig. ^ and compared 
with the lowest order fragmentation function obtained in . A large difference between the 
leading and next-to-leading order quark-to-photon fragmentation functions is observed only 
for z close to 1, indicating the presence of large log(l — z). 

Having described the z distribution in the measured range, we shall now focus only on 
the 'photon' -|- 1 jet rate for z > 0.95 which can be viewed as the isolated photon + 1 
jet rate in the democratic clustering approach 0. This division between isolated and non- 
isolated is motivated both by the measured z distribution (see Fig. ^ and the observation 
that hadronization effects cause the photon jet to have a z value slightly less than unity 0. 
In previous theoretical []T2|, ^ and experimental |^ analyses of isolated photon -|- n jet 



rates, the candidate photon was isolated from the other hadrons using a geometrical cone, 
before these hadrons were clustered into jets by the jet recombination algorithm. The photon 
was said to be 'isolated' if it was accompanied only by a small (but non-zero) amount of 
hadronic energy inside the cone. As a result of these analyses it was found that large negative 
next-to-leading order corrections were needed to obtain a reasonable agreement between 
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Figure 6: The inclusive photon energy distribution normahzed to the hadronic cross section 
as measured by the OPAL Collaboration compared with the 0{a) and 0{aas) order calcu- 
lations including the appropriate quark-to-photon fragmentation functions determined using 
the ALEPH photon + 1 jet data. 



theory and experiment 21, 23 1. Interestingly, a lowest order calculation of the photon 



+ 1 jet rate for the cone algorithm using the process independent photon fragmentation 
function measured by ALEPH and a suitable choice of z to define the isolation yields 
good agreement with the data 0. 

Using the results of our calculation and the fitted quark-to-photon fragmentation func- 
tion, we have determined the isolated photon + 1 jet rate in the democratic approach up to 
0{aas). The result of this calculation compared with data from ALEPH |^ and the leading 
order calculation |T0| is shown in Fig. |[ It can clearly be seen that inclusion of the next-to- 



leading order corrections improves the agreement between data and theory over the whole 
range of ?/cut- It is also apparent that the next-to- leading order corrections to the isolated 
photon + 1 jet rate obtained in this democratic clustering approach are of reasonable size 
indicating a better perturbative stability of this isolated photon definition. 

As a final illustration of the generality of our results, we turn to the inclusive photon 



distribution recently measured by the OPAL Collaboration p5[. Here, jets are not defined 
and, rather than the distribution of electromagnetic energy within the cluster, the energy 
distribution of the photon is considered. OPAL are able to identify photons with energies 
with as little as 10 GeV, and have compared their results with the model estimates of 0, §]. 
They find reasonable agreement with factorization scales fip ^ Mz, although the data is too 
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poor to discriminate between the models. Fig. ^ shows our (scale independent) predictions 
for the inclusive photon energy distribution at both leading and next-to-leading order. We 
see good agreement with the data, even though the phase space relevant for the OPAL data 
far exceeds that used to determine the fragmentation functions from the ALEPH photon + 
1 jet data. 

In summary, we have presented the results of a complete calculation |r3| of the 'photon' 
+ 1 jet rate at 0{aas). Although only next-to-leading order in perturbation theory, this 
calculation contains several ingredients appropriate to the calculation of jet observables at 
next-to-next-to-leading order. In particular, we needed to generalize the phase space slicing 



method of fl^, to more than one theoretically unresolved particle and have analytically 
calculated the following double unresolved configurations : the triple collinear factor, the 
soft /collinear factor and the double single collinear factor. The 'photon' + 1 jet rate was then 
evaluated for a democratic clustering algorithm with a Monte Carlo program using the hybrid 



subtraction method of |T8[ . The results of our calculation, when compared to the data on 
the 'photon' -|- 1 jet rate obtained by ALEPH, enabled a first determination of the process 
independent quark-to-photon fragmentation function at 0{aas) in a fixed order approach. 
As a first application, we have used this function to calculate the 'isolated' photon -|- 1 jet 
rate in a democratic clustering approach at next-to-leading order. The inclusion of the QCD 
corrections improves the agreement between theoretical prediction and experimental data. 
Moreover, it was shown that these corrections are moderate, demonstrating the perturbative 
stability of this particular isolated photon definition. Finally, we have computed the inclusive 
photon energy distribution and found good agreement with the recent OPAL data. 
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